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Abstract 



A new recursive procedure to compute the Zassenhaus formula up 
to high order is presented, providing each exponent in the factorization 
directly as a linear combination of independent commutators and thus 
containing the minimum number of terms. The recursion can be easily 
implemented in a symbolic algebra package and requires much less com- 
putational {effort, both in time and memory resources, than previous algo- 
rithms. In addition, by bounding appropriately each term in the recursion, 
it is possible to get a larger convergence domain of the Zassenhaus formula 
when it is formulated in a Banach algebra. 
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1 Introduction 

Products of exponentials of non-commuting indeterminate variables are of fun- 
damental importance in physics and mathematics. As is well known, the Baker- 
Campbell-Hausdorff theorem states that e"^e^ = e^, with 



Here Zm{X,Y) is a homogeneous Lie polynomial in the non-commuting vari- 
ables X and Y. In other words, is a linear combination (with rational coeffi- 
cients) of commutators of the form [Vi, [V2, . . . , [Kn-i, Kn] • • .]] with Vi G {X, Y} 

*Email: Fernando.Casas@uji.es 
^Email: Ander . Murua@ehu . es 



00 





m=2 



for 1 < i < m. We recall that [X, Y] = XY — YX. The first terms in the series 



(1.1) read explicitly 



Z2 = l[X,Y], Z, = ^[X,[X,Y]]-^[Y,[X,Y]], = ^[X,[Y,[Y,X]]]. 

The expression e"''' e^ = e^ is then properly called the Baker-Campbell-Haus- 
dorff formula (BCH for short) and plays a fundamental role in many fields of 
mathematics (theory of linear differential equations Uj , Lie groups [2j , numerical 
analysis [3]), theoretical physics (perturbation theory [1], Quantum Mechanics 
[S], Statistical Mechanics [6l[7|, quantum computing [8]), control theory (anal- 
ysis and design of nonlinear control laws, nonlinear filters, stabilization of rigid 
bodies f9]), etc. (see [lOj for a comprehensive treatment of the algebraic aspects 
of the BCH formula). 

Although the BCH theorem establishes the precise algebraic structure of the 



exponent Z in (1.1), it does not provide simple ways to compute explicitly this 
series. As a matter of fact, the problem of effectively computing the BCH series 
up to arbitrary degree has a long history, and different procedures have been 
proposed along the years, starting with the work of Richtmyer and Greenspan 
in 1965 (see [11] for a review). Most of the procedures lead to expressions where 
not all the iterated commutators are linearly independent (due to the Jacobi 
identity and other identities appearing at higher degrees). Equivalently, the 
resulting expressions are not formulated directly in terms of a basis of the free 
Lie algebra C{X, Y) generated by the symbols X and Y. Of course, it is always 
possible to write these expressions in terms of a basis, but this rewriting process 
is time consuming and require considerable memory resources. This in addition 
is made more difficult due to the exponential growth of the number of terms 
with the degree m. Recently, a new efficient algorithm has been proposed which 
allows one to get closed expressions for Z^ up to a very high degree in terms 
of both the classical Hall basis and the Lyndon basis of C{X,Y) |11| . 

In the paper dealing with the expansion bearing his name, Magnus [1] cites 
an unpublished reference by Zassenhaus, reporting that there exists a formula 
which may be called the dual of the (Baker-Campbell-) Hausdorff formula. This 
result can be stated as follows. 

Theorem 1.1 (Zassenhaus Formula). Let C{X,Y) be the free Lie algebra gen- 
erated by X and Y. Then, e^^'^ can be uniquely decomposed as 

oo 
n=2 

where Cn{X,Y) £ C{X,Y) is a homogeneous Lie polynomial in X and Y of 
degree n. 

The existence of such a formula is an immediate consequence of the BCH 
theorem. In fact, it is clear that e~^e^"*"^ = e^"*"^, where D involves Lie 
polynomials of degree > 1. Now e~^e^'''^ = e*^^"*"^, where D involves Lie 
polynomials of degree > 2 and the process is repeated again. Induction allows 
one to get the general result. 
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By comparing with the BCH formula it is possible to obtain the first terms 
of the formula (1.2) as 



C2{X,Y) = -l[X,Y], Cs{X,Y) = ^[Y, [X,Y]] + ^[X, [X,Y]]. 

Although less familiar than the BCH formula, the Zassenhaus formula con- 
stitutes nevertheless a standard tool in several fields, since the disentangling of 
the exponential of the sum of two non-commuting operators into an (in gen- 
eral infinite) product of exponential operators arises for instance in statistical 
mechanics, many-body theories, quantum optics, path integrals, g-analysis in 
quantum groups, etc. [12]. Also in particle accelerators physics, Dragt and 
his collaborators have used the Zassenhaus formula to compute the relevant 
maps both in Taylor and factored product form |13j . Yet in another context, 
very recently Iserles and Kropielnicka [13] have proposed a new family of high- 
order splitting methods for the numerical integration of the time-dependent 
Schrodinger equation based on a symmetric version of the Zassenhaus formula. 

Several systematic computations of the terms C„ for n > 3 in the Zassenhaus 
formula have been carried out in the literature, starting with the work of Wilcox 
|7| , where a recursive procedure is presented that has been subsequently used to 
get explicit expressions up to Cg in terms of nested commutators [12J. On the 
other hand, Volkin [15] proposed a general technique to express a function of 
non-commuting operators as an expansion in successively higher commutators 
of the operators involved. In particular, he was able to get recursive formulae 
up to Cg. By following an idea already suggested by Wilcox in [tj, Suzuki [T6] 
obtained the successive terms Cn{X,Y) in 

by differentiating both sides with respect to A and setting A = after each 
differentiation. In this way 

and in general, for n > 3, 

n! VrfA"^ Va=o 

Finally, Baues |17j gave explicit formulae for the Zassenhaus terms via ho- 
motopy theory and more recently Kurlin |18j obtained a closed expression for 
J2n>2 ^ri in the metabelian case. 

All of these proposals give results for Cn as a linear combination of nested 
commutators. In contrast, Scholz and Weyrauch [19] proposed a recursive pro- 
cedure based on upper triangular matrices that can be easily implemented in 
a symbolic algebra package. In this case, however, the expressions for C„ are 
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not explicitly written down in terms of homogeneous Lie polynomials. More 
recently [20], the same authors have applied a technique related to one previ- 
ously used by Witschel [21] to get Cn up to n = 15 in less than 2 minutes of 
CPU time. Here also the Zassenhaus terms are written as 

Cn= ^ 9wW, (1.5) 
ui,|ui|=?i 

where is a rational coefficient and the sum is taken over all words w with 
length = n in the symbols X and Y , i.e., w = 0102 • • • a„, each being 
X or Y . Of course, one may always apply the Dynkin-Specht-Wever theorem 
[22], and express Cn as 

Cn= - 9w [w], (1-6) 

w,\w\=n 



that is, the individual terms are the same as in the associative series ( 1.5 ) except 
that the word w = aia2 . . . a„ is replaced with the right nested commutator 
[w\ = [ai, [a2, . . . [an_i,an] . . .]] and the coefficient g^j is divided by the word 
length n. In this way Cn is constructed as a linear combination of nested 
commutators of homogeneous degree n, that is, as a linear combination of 
elements of the homogeneous subspace C{X^ Y)n of degree n of the free Lie 



algebra C{X,Y). As a matter of fact, another representation of (1.5) in terms 
of nested commutators is proposed in ^20j which, it is claimed, contains less 
terms than the Dynkin-Specht-Wever representation. In any case, it should 
be stressed that the set of nested commutators [w] for words w of length n in 
either representation is not a basis of the homogeneous subspace C{X,Y)n. 

The purpose of this paper is twofold. First, to present a new recurrence 
that allows one to express the Zassenhaus terms C„ directly as a linear com- 
bination of independent elements of the homogeneous subspace C{X,Y)n- In 
other words, the procedure, which can be easily implemented in a symbolic 
algebra package, gives C„ up to a prescribed degree directly in terms of inde- 
pendent commutators involving n operators X and Y . In this way, no rewriting 
process in a basis of C{X, Y) is necessary, thus saving considerable computing 
time and memory resources. Moreover, we are able to express directly Cn with 
the minimum number of commutators required at each degree n. 

The second aspect we are dealing with concerns the convergence of the 
Zassenhaus formula when it is formulated in a Banach algebra. As far as we 
know, there are only two previous results in the literature. The first one was 



obtained by Suzuki [T^ starting with the recurrence (1.4). Specifically, he was 



able to prove that if |A|(||X|| + \\Y\\) < log 2 - 1/2 the infinite product (1.3) 



converges. Subsequently, Bayen [23] generalized the analysis, showing that the 



product (1.3) converges if |A|(||X|| -|- \\Y\\) < r, where r 0.596705 is a root 



of a certain transcendental equation. In the present work, we obtain sharper 
bounds for the terms of the Zassenhaus formula which show that the product 



(1.2) converges in an enlarged domain. 

A simple but important remark is in order here. In some applications, the 
"left-oriented" Zassenhaus formula 

e^+^ = . . . e'^^(^'^) e^^(^'^) e'^^^^'^) e^ e^ (1.7) 
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is used instead of (1.2). A simple observation shows that the exponents Ci and 



Ci are related through 

Ci{X,Y) = {-iy+^Ci{X,Y), i>2 



and so we may restrict ourselves to analyzing the "right-oriented" formula ( 1.2 ). 



2 Constructing the Zassenhaus terms 
2.1 A new recurrence 

To derive our recursive procedure, it is convenient to introduce a parameter A 

^HX+Y) ^ ^XX ^XY ^X^C, gA3C3 gA4C4 . . . (2.I) 



as in (1.3) 



so that the original Zassenhaus formula (1.2) is recovered when A = 1. More- 
over, we consider the compositions 



i?i(A) = e-'^e 



and for each n >2, 
i?„(A)=e-^"^" 
It is then clear that 



g-A^Ca g-AF g-AX ^X{X+Y) ^ g-A"C„ 



(2.2) 



i?„_i(A). (2.3) 



(2.4) 



Finally, we introduce 

Fn{X) 



Rn{X)]Rn{Xr\ 



n > 1. 



(2.5) 



To determine the series Fn{X) we proceed as follows. On the one hand, a simple 
ition st 

FniX) 



calculation starting from (2.3), leads for n > 2 to 

d 



-nC„A"-^ +e- 



dX 



i?„_i(A))i?„_i(A)-ie^"^" 



-nC„A"-^+e-^"^"F„.i(A)e^"'^" 

-naA"-^+e-^""<^^"F„_i(A) 

-^"'^c^{Fn-i{X)-nCny'''), 



(2.6) 



where we have used the well known formula 



n>0 



with 



adA^ = [A, B], ad^B = [A, ad^"^S], ad^B = B. 
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On the other hand, differentiating expression (2.4) with respect to A and 
taking into account (2.5) we arrive at 

oo 

F„(A) = (n + l)C„+iA"+ ^ jA^-ie""^''^'^^"+i---e"'"'^'^^.-C,. (2.7) 

j=n+2 



In other words, 



F„(A) = (n + l)C„+iA" + G,+i(A), n > 1, 



where G„+i(0) = Gi'|i(0) = ■ 
n > 1, 

Fn+l{\) 
Cn+l 



G^"i^^(0) = 0. In consequence, we have, for 



e-A"+^adc„^, G„+i(A), 
A*^ 



G„+i(A) = F„(A)--F(")(0). 

n! 



(2.8) 
(2.9) 

(2.10) 



Expressions (2.8)-(2.10) ahow one to compute recursively the Zassenhaus terms 
Cn starting from -Fi(A). The sequence is 

Fn{X) Cn+l Gn+i{\) F„+i(A) n > 1. 

Let us analyze more in detail this procedure, with the goal of providing an 
algorithm well adapted from a computational point of view. 
For n = 1, and taking into account (2.2), we get 

Fi(A) = 



Y - e-^^Xe^^ + e-^^e-"^(X + y)e"^e^^ 



-Aady ( —\a.dx 



that is, 



i=0 j=l 



(2.11) 



or equivalently 

oo 



with = ^ 



^ ^^'^ -ad^-^ad^;^y. (2.12) 



k=l 



In general, from (2.6) a straightforward calculation shows that for n > 2, 

oo [k/n]-l 



Fn{X) = ^^fn,kX^, with fn,k= ^ —&<5^(j^fn-l,k-nj-, 



k > n. 



k=n 



3=0 



(2.13) 
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Here [k/n] denotes the integer part of k/n. Moreover, a closer examination of 



(2.7) reveals that 



(n + l)C„+iA" + (n + 2)e^"^'"'^^"+i C„+2A"+i + 

(n + l)Cn+iy' + (n + 2)C„+2A"+^ + • • • 

+(2n + 2)C2„+2A2"+i + A2«+2[c„+i, C7„+2] + • • • 

2n+2 

^ A:CfcA'=-i + A2"+2^„(A), 

fc=n+l 



(2.14) 



where Hn{\) involves commutators of Cj, j > + 1. Notice that the terms 
Cn+ii . . . , C2n+2 of the Zasscnhaus formula can be then directly obtained from 



F„(A). In particular, one directly gets from (2.11) 

ra-l 

~\{n-j)\ 

Explicitly, 



1 """^ (-^\ 

V-) ^radtad^-^y, for n = 1, 2, 3. (2.15) 



C2 
C4 



1 



-i([y, [y, [x,y]]] + [y, [x, [x,y]]]) - ^[x, [x, [x,y]]]. 



1 



Taking into account (|2.13|) and (|2.14|) we have in general 

n > 5, 



Cn 



1 



+1 



n 



1 f[n/2], 



(2.16) 



where the expressions of /„ ^ are given recursively by (2.13). 



In summary, the algorithm we propose for computing the Zassenhaus terms 
is the following: 



Define by eq. ^A2) 

Cn = (l/n)/i,„_i, n = 2,3,4 

Define fn,k n>2, k >n by eq. (2.13) 

Cn = (l/?^)/[(n-l)/2],n-l n > 5. 



(2.17) 



This constitutes a new recursive way for obtaining directly the term C„ as a 
homogeneous Lie polynomial in X, y of arbitrarily large degree n which can 
be easily implemented with a symbolic algebra package. 

2.2 Linear independence 



Algorithm (2.17), or equivalently the procedure given by the identities (2.8) 



(2.10), provides expressions for C„ that, by construction, involve only indepen- 



dent commutators. In other words, they cannot be simplified further by using 
the Jacobi identity and the antisymmetry property of the commutator. 
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In order to prove this assertion, it is convenient to get a more explicit 
expression of Fn{X) and Cn+i from (2.8)-(2.10). To this end, consider for 
n > 1 the sets Jn and In of (n + l)-tuples of non- negative integers recursively 
defined as follows: 



Ji = {(^o,ii) eN^ : to > 1}, 

= {{ioik, ■ ■ ■ ,in) ^ Jn ■ io + H + H \-nir. 

Jn+1 = {Jn\In) X N. 



n}. 



The set X„ can be directly defined as the set of [n + l)-tuples of non- negative 
integers satisfying that iq + ii + 2i2 + ■ ■ ■ + nin = n and 



j + 1 <iQ + ii+2i2 h i i,- for j = 0, 



1. 



(2.18) 



Thus, in particular, Zi = {(1,0)}, Z2 = {(1, 1, 0), (2, 0, 0)}, etc. Observe that, 
by construction, each (io, ii, in) £ 2r„ is such that im = if m > n/2. 
From (2.8)-( [2ri0| ), one can then prove by induction on n that, for n > 1, 

(_X)*0~l hin^io+«l+2i2H h' 



n+1 



E 

(i0,jl,---,in)Gj'n 
1 

n + 1 



E 



ioHi! • • - in! 

_]^yoH \-in 



-ad 



(i0,il,...,jn)£ln 



-ad*;? 



••adg^ad^ad^F. 



cLCi^^ clCiy^ ctdjY" ^5 

(2.19) 



In fact, this is clearly true for n = 1 (equations (2.11) and (2.15), respectively), 
whereas successive application of (2.8)-(2.10) leads to the general result. 

Now, repeated application of Lazard elimination principle |24j . together 
with Xi = {(1,0)}, {C2} = {-l[X,Y]} = {-iad^^ad^y : {io,h) G XJ, 
shows that, vector space, 

C{X,Y) = span({X}) e £({ad^y : j > 0}) 

= span({X}) e span({y}) £({ad*yad*^y : i > 0, j > 1}) 

= span({X,y})e£({ad^ad*°y : {io,k) e Ji}) 

= span({X,y})e£({C2}U{ad^ad^y : {io,ii) e Ji\Ii}) 

= span({X,y,C2})©£({adg^ad*^ad^y : (zq, ii, ^2) G ^^2}). 

More generally, application of Lazard elimination together with ( |2.19| ) gives 

£(X,y) C span({X,y,C2,...,C„}) 

e/:({ad^^---adg^ad*^ad^y : (io, • • • , in) G Jn}) 
C span({X,y,C2,...,C„}) 

© £({Cn+i} U {ad^^ • • • adg^ad^ad^y : (io, . . . , in) G Jn\In}) 
C span({X,y,C2,...,C„+i}) 



}C{{ad 



Cn + 1 



adg^adyad^y : (io, . . . , i„+i) G Jn+i})- 



In consequence, the terms {ad^ • • • ad^^adyad^y : (io, ii, • • • , im) G Jn} 
are linearly independent in the free Lie algebra C{X,Y) and the same is true 
for the representation (2.19) of the Zassenhaus terms. 



8 



2.3 Computational aspects 



We have implemented the recursive procedure (2.17) in Mathematical^ as the 
following algorithm. 

Clear[Cmt, ad, ff, cc] ; 
$RecursionLimit= 1024; 
Cmt [a_ , a_] : = ; 

Cmt [a 

Cmt [a 

Cmt [a 



0, b ] := 0; 

c_ + d_, b ] := Cmt[a, c, b] + Cmt[a, d, b] ; 

n_ c_Cmt, b ] := n Cmt [a, c, b] ; 

n_ X, b ]:= n Cmt [a, X, b] ; 

n_ Y, b ]:= n Cmt [a, Y, b] ; 



Cmt [a_. 
Cmt [a_. 

Cmt /: Format [Cmt [a_ , b_]]:= SequenceForm[" [" , a, b, "] "] ; 

ad[a_, 0, b_] := b; 

ad[a_, j_Integer, b_] := Cmt[a, ad[a, j-1, b]] ; 
ff [1, k_] := ff [1, k] = 

Sum[((-l)-k/(j ! (k-j)!)) ad[Y, k- j , ad[X, j, Y]], { j , 1, k}] ; 
cc[2] = (1/2) ff[l, 1]; 
ff [p_, k_] := ff [p, k] = 

Sum[((-l)-j/j !) ad[cc[p], j, ff[p-l, k - p j]], {j , 0, 
IntegerPart [k/p] - l}] ; 
cc [p_Integer] : = cc [p] = 

Expand[(l/p) ff [IntegerPart [(p-l)/2] , p-1]]; 

The object Cmt[xi, X2, . . . , x„_i, Xn] refers to the nested commutator 
[xi, [j;2, • • • [2;n-i, a^n] • • • ]]• It has attached just the linearity property (there 
is no need to attach to it the antisymmetry property and the Jacobi iden- 
tity). The symbol ad represents the adjoint operator and its powers ad^6, 



whereas ff [l,k], ff [p,k] and cc[p] correspond to expressions (2.12), (2.13) 



and (2.16), respectively. Proceeding in this way, we have obtained the expres- 
sions of Cn up to n = 20 with a reasonable computational time and memory 
requirements. Thus, for instance, constructing the terms up to degree n = 20 
with a personal computer (2.4 GHz Intel Core 2 Duo processor with 2 GBytes 
of RAM) takes less than 20 seconds of CPU time and 35 MBytes of memory. 
The expression for C20 has 48528 terms, all of them independent. The resulting 
expressions up to Cg are identical to those expressed in the classical Hall basis. 

In Table [1] we collect the CPU time (in seconds) and memory (in MBytes) 
needed to construct the terms C2, C3, . . . , up to a given value of n both with 



the recurrence ( |2.17 ) (New) and the implementation provided in [20j using a 



variant of the so-called comparison method previously introduced in [21J ( W-S). 
Notice that with this method, which is the most efficient of all the procedures 
analyzed in [20], the Zassenhaus exponents Cj are expressed as linear combina- 
tions of words of length j and not directly in terms of independent commutators 
(although this is always possible by applying Dynkin-Specht-Wever theorem or 
Theorem 2 in [20], and then simplifying the resulting expressions by taking into 
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account the Jacobi identity and the antisymmetry property of the commuta- 
tor). For comparison, Cig has 54146 terms when expressed as combinations of 
words, but only 3711 terms with the new formulation. This translates directly 
into the memory requirements of both algorithms, as is evident from the results 
collected in the table. 



n 


CPU time (seconds) 


Memory (MBytes) 




W-S 


New 


W-S 


New 


14 


29.18 


0.14 


122.90 


0.88 


16 


203.85 


0.59 


764.32 


4.09 


18 




3.01 




11.12 


20 




19.18 




35.27 



Table 1: CPU time and memory required for the computation of the Zassenhaus terms 
C2 , C3 , . . . , C„ up to the given value of n using the algorithm presented in [5D] ( W-S) 
and recurrence (2.17) {New). 



3 Convergence of the Zassenhaus formula 

Suppose now that X and Y are defined in a Banach algebra A, that is to 
say, an algebra that is also a complete normed linear space whose norm is 
submultiplicative , 

||xy|| < ||x|| ||y|| (3.1) 

for any two elements of A. Notice that for the commutator one has || [X, Y]\\ < 
2 ||X|| \\Y\\. Then it makes sense to analyze the convergence of the Zassenhaus 



formula (1.2). 

As stated in the introduction, we are aware of only two previous results 
establishing sufficient conditions for convergence of the form ||X|| + ||y|| < r 
with a given r > 0. Specifically, Suzuki [16] obtained Vg = log 2 — 2 ~ 0.1931, 
whereas Bayen [23] proved that the domain of convergence can be enlarged up 



to a value of given by the unique positive solution of the equation 

zMl + 2y — - — dwj = 4(21og2 - 1). 

A numerical computation shows that rb = 0.59670569.... Thus for ||X|| + 
||y|| < rf, one has 

lim e^e^e^^---e^"=e^+^. (3.2) 



In the following, we use recursion (2.8)-(2.10) to show that (3.2) holds 
indeed for {x,y) = (||Ar||, \\Y\\) G in a domain that is larger than {(x,y) G 
: 0<x + y < n}. 



Clearly, (|3^ holds if 

lim WRnim = 1, (3.3) 
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where Rn{^) is given by (2.4), and thus is the solution of the initial value 
problem 

^Rn{X) = Fn{X)Rn{X), Rn{0) = I. (3.4) 

It is well known that, if ||F„(A) ||(iA < oo, then there exists a unique solution 
Rn{X) of Q for < A < 1, and that ||i?„(l)|| < exp(/(J ||F„(A)||(iA). In con- 
sequence, convergence (3.3) will be guaranteed whenever = (||^||, ||^||) £ 
is such that 



lim 

n— >oo 



|F„(A)||dA = 0. 



From (2.19) we have that ||C„+i|| < 6n+i, where 62 = xy and for n > 2, 

2*0 H hin 



n + 



1 E 



{io,jl,.--i*n)Glr 

Similarly, ||F„(A)|| < /n(A), where 

(2A) 



ioliil ■■■ij. 



il=0 io=l 



(3.5) 



and for n > 2, 

/n(A) = 

Note that this implies 



E 



2*0 H hin_>^«0+n+2i2H \-r. 



„i 00 

/ /n(A)(iA< V<5fc, 
•^0 k=n 



SO that (3.3) is ensured if the series ^^=2^k converges. Let us analyze each 
term of this series by mimicking the recursive procedure given by (2.17). From 
(2.12) (or alternatively from (3.5)) and (2.13), we get 



[fc/n]-l 



~^^n^ri—l,k—nj 



whence 



\Cn\\ <Sn = -C^[(n-l)/2],n 



n 



-1) 



n > 3. 



A sufficient condition for convergence is obtained by imposing 

lim %ti < 1. 

n— ^00 On 



(3.6) 



(3.7) 



(3.8) 
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At this point it is worth remarking that, although not reflected by the notation, 
both dn,k and (5„ depend on = (||^||, \\Y 

fact a constraint on the convergence domain {x,y) G 



so condition (3.8) imphes in 
^ of the Zassenhaus 



formula. In Figure [T| we show graphically the (numerically computed) domain 
D of such points {x,y). This has been obtained by computing for each point 
the coefficients dn,k and 5„ up to n = 1000 (in fact, considering a smaller 
value of n the figure does no change significantly). We have also included 
for comparison the previous results x + y < 0.1931 and x + y < 0.5967 of 
Suzuki and Bayen, respectively. Observe that the new convergence domain is 
considerably larger. In particular, it contains the region x + y < 1.054, but it is 
not restricted to that. For instance, the convergence domain contains the sets 
{(x, 2.9216) : x < 0.00292} and {(2.893, y) : y < 0.0145}, and also the points 
(x, 0) and (y, 0) with arbitrarily large value of x or y. 



3.0 



2.5 



2.0 




Figure 1 : Upper boundary of convergence domain for the Zassenhaus formula obtained 



with the procedure (13J|)-(13J|). Previous results < ||X|| + < 0.1931 and < 
||X|| + ||y|| < 0.5967 are also depicted for comparison. The new domain contains the 
region \\X\\ + \\Y\\ < 1.054. 

In summary, we have presented a new recursive procedure that not only 
allows us to get the expressions of the Zassenhaus exponents Cn directly in terms 
of independent commutators in an efficient way and can be easily implemented 
in any symbolic algebra package, but also shows, by bounding appropriately 
each term in the recursion, that the convergence domain of the Zassenhaus 
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formula is considerably larger than the domain guaranteed by previously known 
results. 
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